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Abstract. A generally relativistic theory of thermodynamics is developed, based on four 
main physical principles: heat is a local form of energy, therefore described by a thermal 
energy tensor; conservation of mass, equivalent to conservation of heat, or the local first 
law; entropy is a local current; and non-destruction of entropy, or the local second law. A 
fluid is defined by the thermostatic energy tensor being isotropic. The entropy current is 
related to the other fields by certain equations, including a generalised Gibbs equation for 
the thermostatic entropy, followed by linear and quadratic terms in the dissipative (thermal 
minus thermostatic) energy tensor. Then the second law suggests certain equations for 
the dissipative energy tensor, generalising the Israel-Stewart dissipative relations, which 
describe heat conduction and viscosity including relativistic effects and relaxation effects. 
In the thermostatic case, the perfect-fluid model is recovered. In the linear approximation 
for entropy, the Eckart theory is recovered. In the quadratic approximation for entropy, the 
theory is similar to that of Israel & Stewart, but involving neither state-space differentials, 
nor a non-equilibrium Gibbs equation, nor non-material frames. Also, unlike conventional 
thermodynamics, the thermal energy density is not assumed to be purely thermostatic, 
though this is derived in the linear approximation. Otherwise, the theory reduces in 
the non-relativistic limit to the extended thermodynamics of irreversible processes due to 
Miiller. The dissipative energy density seems to be a new thermodynamical field, but also 
exists in relativistic kinetic theory of gases. 
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I. Introduction 

Thermodynamics is perhaps best known to many relativists by its analogies with black- 
hole dynamics, made concrete by the famous result of Hawking that stationary black holes 
radiate quantum fields with a thermal spectrum. Recently, another link has been obtained: 
a unified first law which contains first laws of both black-hole dynamics and relativistic 
thermodynamics, as described in an accompanying paper [1]. This result requires lit- 
tle thermodynamics other than the general form of the energy tensor given long ago by 
Eckart [2]. However, the literature does not seem to contain a fully satisfactory general 
theory of relativistic thermodynamics, in the classical sense of a general macroscopic theory 
independent of (but consistent with) microscopic theories. This remarkable situation has 
prompted the development of the theory described in this paper. Such a theory naturally 
has widespread applications in astrophysics and cosmology. For simplicity, only a single 
material will be considered here, ignoring mixtures, phase changes, chemical reactions and 
external forces. 

Perhaps the greatest progress in relativistic thermodynamics was made by Eckart 
[2], who: identified the basic fields, formulated in this paper as the material current, the 
thermal energy tensor, the entropy current and the temperature; gave the basic conser- 
vation equations, of mass and energy-momentum; related the entropy density for fluids 
by a relativistic Gibbs equation; and gave a relativistic second law, showing that it could 
be satisfied by simple dissipative relations which generalise the Fourier and Navier-Stokes 
equations for heat conduction and viscosity. Thus the thermodynamic system is described 
by a finite number of macroscopic fields which may be measured by standard techniques. 

Although Eckart was concerned with special relativity, the theory generalises imme- 
diately to general relativity. The only unacceptable part of the Eckart theory is that 
the dissipative relations entail non-causal propagation of disturbances in temperature and 
velocity. The problem of infinite propagation speeds existed even in non-relativistic ther- 
modynamics for longer still, until Miiller [3] showed how to resolve it by allowing the 
entropy current to contain terms which are quadratic (rather than just linear) in the dissi- 
pative fields, namely thermal flux and viscous stress. Physically, this means assuming that 
the dissipative fields are small enough to justify the approximation: thermostatic, linear 
or quadratic. This incidentally means that entropy is generally not a state function, as 
claimed in classical thermodynamics, though this is retained in the linear approximation. 

One might expect to obtain an acceptable theory of relativistic thermodynamics sim- 
ply by applying a Miiller extension to the Eckart theory. This was attempted by Israel 
& Stewart [4], who obtained modified dissipative relations which are consistent with the 
Miiller theory. However, they gave a type of Gibbs equation which differs from that of 
Eckart and does not reduce appropriately to the Gibbs equation of non-relativistic local 
thermodynamics, as described for example by de Groot & Mazur [5]. Also, Israel & Stew- 
art gave a theory intended to hold not only in the material frame of Eckart, but also in the 
frame of Landau & Lifshitz [6] . This seems to involve some further approximation which 
is unclear, at least to this author. In any case, this paper will use the material frame only, 
for reasons explained below. 

The theory described in this article was developed from first principles, but includes 
Muller-type quadratic dissipative effects, without which the Eckart theory is recovered. 
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The theory is based on simple physical principles which are consistent with local non- 
relativistic thermodynamics, including local first and second laws. These principles serve 
to define what is meant by thermodynamic matter, which is intended to describe real fluids 
and solids in circumstances where a macroscopic description is adequate. The most basic 
principle is that heat is a local form of energy. In relativity, heat is therefore described by 
an energy tensor, called the thermal energy tensor, whose various components are thermal 
energy density, thermal flux and thermal stress. Thermal flux is what is normally called 
heat flux only in the material frame, thereby fixing this frame as the physically natural 
one. Thermal stress is what is normally called material stress. Thermal energy is what is 
normally called internal energy, but with fixed zero. 

This principle entails one further generalisation: the thermal energy density need 
not be purely thermostatic, allowing a further dissipative field, namely dissipative energy 
density. Classical thermodynamics excludes this by the claim that internal energy is a state 
function. However, this is inconsistent with relativistic kinetic theory of gases [7], where 
thermal energy density is defined in terms of the molecular distribution and is generally not 
purely thermostatic. In fact, dissipative energy density emerged unrecognised even in non- 
relativistic kinetic theory as essentially the 14th moment of the 14-moment approximation 
of Grad [8]. It turns out that relativistic causality requires dissipative energy density 
to vanish in the linear approximation for entropy, thereby recovering the Eckart theory. 
However, this need not hold in the quadratic approximation, thereby generalising the 
Israel- Stewart dissipative relations. 

This article is intended for relativists who may have no previous experience of thermo- 
dynamics. (Thermodynamicists with a nodding acquaintance with relativity should also 
find it accessible). Therefore brief reviews are given of Miiller's extended thermodynam- 
ics, of the previously standard local (or field) thermodynamics, and of relevant classical 
thermodynamics. This seems necessary because most textbooks on thermodynamics, even 
today, insist on applying thermostatic concepts in ways which are quite inapplicable in gen- 
uine thermodynamics, even in the formulation of the first and second laws. Such confusions 
often survive in work on relativistic thermodynamics. 

The article is organised as follows. Sections II- VI respectively review classical ther- 
modynamics, non-relativistic hydrodynamics, non-relativistic thermodynamics, extended 
thermodynamics and basic general relativity. Section VII introduces the relativistic con- 
cept of heat and the local first law, or conservation of heat, equivalent to conservation of 
mass. Section VIII checks the Newtonian limit. Section IX divides the thermal energy ten- 
sor into thermostatic and dissipative parts and defines a fluid by the thermostatic energy 
tensor being isotropic. Section X introduces the relativistic concept of entropy and the 
local second law, and describes the thermostatic case using a generalised Gibbs equation. 
Section XI describes the linear approximation for entropy, showing that it reduces to a mi- 
nor generalisation of the Eckart theory. Section XII describes the quadratic approximation 
for entropy, obtaining dissipative relations which generalise the Israel-Stewart relations (in 
the material frame). Section XIII concerns integral quantities and corresponding laws. 
Section XIV concludes. A point of style is that equations in the text are sometimes used 
for temporary arguments, whereas displayed equations are intended to carry more weight. 
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II. Classical thermodynamics 

In classical thermodynamics, the basic quantities are temperature heat supply Q, work 
W, internal energy H and entropy S. The classical first law is said to be dH = 5Q + 5W 
and the second law $dS > 5Q, where d is said to be a differential and 5 not, neither 
being actually defined except in thermostatics, where the second law becomes an equality, 
or in the limit where they become difference operators between two equilibrium states 
separated in time, corresponding to the classical experimental situation. Clearly this is 
nonsense, as has been emphasised eloquently by Truesdell [9]. The resolution is simply 
that these operators are really derivatives with respect to time. In fluid mechanics the 
relevant derivative is that along the fluid flow, called the material or comoving derivative, 
henceforth denoted by a dot. Similarly, in solid mechanics the relevant derivative is that 
in the centre-of-momentum frame, which may also be called the material derivative. Then 
the first law is 

H = Q + W (2.1) 

where Q and W should be taken as the basic quantities, Q and W themselves not being 
uniquely defined. For instance, for an inviscid fluid, the work is given by W — —pV, where 
V is the volume and p the pressure of the fluid, so that the first law reads 

H = Q-pV. (2.2) 

Similarly, the second law is 

•&S > Q. (2.3) 

Both laws now make sense, though it turns out that they respectively require the pres- 
sure and temperature to be uniform, again reflecting the classical experimental situation. 
Here and henceforth, uniform means spatially constant. In general, the laws need to be 
formulated locally, as described subsequently. 

The classical zeroth law states that temperature $ is constant in thermal equilibrium. 
The classical third law, expressed in its loosest form, is that absolute zero temperature 
cannot be physically attained. Thus ■d > in practice. 

The entropy may be divided into entropy supply S , given by 

$S = Q (2.4) 

and entropy production S — S a - Then the second law may be written 

S > S (2.5) 

which expresses entropy production. In words, S is the entropy of the system, where system 
means a comoving volume of material, and S is the entropy supplied to the system. Thus 
the second law implies that the total entropy of the universe cannot decrease. 
Equality in the second law, 

S = S (2.6) 
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holds in thermostatics, often called equilibrium thermodynamics or reversible thermody- 
namics. In the thermostatic case, the first and second laws for an inviscid fluid imply 

&S = H + pV (2.7) 

which is the Gibbs equation. 

It was found experimentally that the thermostatic state of a gas could be described 
by just two quantities, such as density p and temperature. For uniform density, one may 
use the material volume V instead of density, since mass M = V p is conserved: 

M = 0. (2.8) 

In particular, H and p are said to be state functions, meaning functions of (V, ■#). If S is also 
a state function, then the Gibbs equation can be rewritten in terms of the differential d in 
state space as fidS = dH+pdV, as Gibbs originally did [10]. This follows because the state 
functions / are now uniform by assumption, so that df = fdt. However, in the non-uniform 
case, the last relation does not hold and the two types of Gibbs equation cannot both be 
generally correct. It turns out that the material Gibbs equation (2.7), localised, forms 
part of the thermodynamic field equations [5], but that the state-space Gibbs equation is 
unnecessary.* This point deserved emphasis because the Gibbs equation is often given in 
state-space form in introductory textbooks on thermodynamics, or on mechanics including 
thermodynamics, as well as in much work on relativistic thermodynamics. Indeed, this 
misunderstanding seems to be the root of the meaningless first and second laws mentioned 
above. The thermodynamic theory developed in this article will involve the material 
derivative rather than state-space differentials. 

III. Non-relativistic hydrodynamics 

In non-relativistic hydrodynamics the basic quantities are the velocity vector v, density p 
and stress tensor r of the fluid. The stress tensor is divided into thermostatic pressure p 
and viscous stress a by 

T = a +ph (3.1) 

where h is the flat spatial metric. Some authors use the opposite sign convention for stress, 
but the above convention is standard in relativity. The viscous stress a may be further 
divided into trace and traceless parts 

w = \h \ a (3.2a) 
q = a - \{h : a)h (3.26) 

* For instance, in a static, spherically symmetric case with radius r, one has df = f'dr 
and the state-space Gibbs equation implies $5" = H'+pV . For liquid helium II at absolute 
zero, this implies an unphysical equation of state H = —pV, whereas the material Gibbs 
equation (2.7) would relate incompressibility, V = 0, to conserved internal energy, H = 0. 
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where the colon denotes double symmetric contraction. Then w is the viscous pressure (or 
bulk viscous stress) and q is the shear (or deviatoric) stress, satisfying 

a = q + zuh. (3-3) 

The inviscid case is defined by a = 0. 

Denoting the spatial volume form by * and a spatial region by E, the volume is 
recovered as 

V= f *1 (3.4) 



s 



and the mass as 



Note the kinematic relation 



M 



[ *P- (3-5) 



* = *D ■ v (3-6) 



where D is the spatial gradient and the centred dot denotes contraction. Then conservation 
of mass (2.8) may be written in a local form as the continuity equation 

p + pD ■ v = 0. (3.7) 

Similarly, conservation of momentum, or Newton's first law for the fluid, yields the Euler- 
Cauchy equation 

p(v + D&) = —D ■ t (3.8) 

where $ is the Newtonian gravitational potential. Conservation of mass and momentum 
are the basic hydrodynamic field equations. 

IV. Non-relativistic thermodynamics 

The classical theory of thermodynamics must be reformulated locally to be compatible 
with hydrodynamics. This can be done by introducing local quantities which replace the 
classical integral quantities. Firstly there are the internal energy density e and entropy 
density s, in terms of which 



*e (4.1a) 

S= I *s. (4.16) 
Jt, 

Actually it is more traditional to use the specific internal energy e* = e/p and the specific 
entropy s* = s/p, due to frequent occurrences of combinations of the form 

pf* = f + fD-v (4.2) 

where /* = f/p for a function /; the identity follows from conservation of mass (3.7). 
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Heat supply is replaced locally by a heat flux vector g, in terms of which 



Q 



-I .q = -*Dq (4.3) 

JdT, JE 



where the second expression follows from the Gauss divergence theorem. Similarly, defining 
the entropy flux vector <p = q/tf allows the entropy supply to be recovered as 

So = -<f -<p= - f *D-<p. (4.4) 

This relates to Q as in the classical definition (2.4) if the temperature is uniform, D-d = 0. 
The local first law is 

e + eD-v = -D-q-r : (D ® v) (4.5) 

where <g> denotes the symmetric tensor product. This local first law integrates to the 
classical form (2.2) for an inviscid fluid with uniform pressure, Dp = 0. Similarly, the local 
second law is 

s + sD-v + D-ip>0 (4.6) 

which integrates to the classical second law (2.5). Lastly, the local form of the Gibbs 
equation is 

tfs* = e* -pp/p 2 - (4.7) 

For uniform pressure and temperature, this integrates to the classical Gibbs equation (2.7). 

It was found experimentally that the thermostatic state of a gas could be described 
by the density and temperature only. Often this is generalised by giving two equations 
relating (p, e, p) which specify a two-dimensional subspace. However, this article will 
maintain (p, $) as the state variables for definiteness. In particular, p and e* are functions 
of (p,$), these relations being called, respectively, thermal and caloric equations of state. 
The simplest equations of state are those of an ideal gas: 

p = Ropti (4.8a) 
e* = c tf (4.86) 

where Rq is the gas constant and cq is the specific heat capacity at constant density, 
usually called constant volume. The local Gibbs equation then integrates to give the 
specific entropy, up to a time-independent function, as 

s* = c log ?? - #0 log p. (4.9) 

In the inviscid case, the first law yields 

-D q = c p$ - Rotip. (4.10) 
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Comparing with the classical definition of heat capacity Q/$, implicitly assuming uniform 
temperature, it can be seen that the specific heat capacities at constant density and pres- 
sure are cq and c$ + Rq respectively. Experimentally, Rq and cq are nearly constant at 
low pressure, for a wide range of density and temperature [11]. Moreover, it is found that 
there is a universal gas constant Rotuq independent of the type of gas, where mo is the 
molecular mass, and that monatomic gases satisfy 

c = §#o- (4.11) 

This relation and the ideal gas laws (4.8) can be derived in the kinetic theory of gases [5], 
with RqUIo being the Boltzmann constant. 

Completing the system of equations requires further equations for (q, a) which are 
consistent with the second law. This is a tight restriction, as may be seen by eliminating 
the internal energy between the first law (4.5) and Gibbs equation (4.7), yielding 

q-Dti <t:(D®v) 
s + sD.v + D^ = -^- 2 y - i. (4.12) 

The simplest way to ensure compliance with the second law (4.6) is for the right-hand side 
to be a sum of squares, which leads to 

q = -K Dd (4.13a) 

w = -X D-v (4.136) 

q = -2hq{D ® v - \{D ■ v)h) (4.13c) 

with 

k > A > /U > 0. (4.14) 
If these coefficients are non-zero, the entropy production is given explicitly by 

2 

s + sD . v + D .^ = l^ + ™ + llS_>0 (4.15) 

which is manifestly non-negative. This makes it clear that (q, a) cause entropy produc- 
tion, i.e. they are thermally dissipative in nature. In this article, (q, a) will be called the 
dissipative fields and the equations (4.13) determining them the dissipative relations. The 
dissipative relations and the equations of state are usually described collectively as the 
constitutive relations. 

One may recognise the dissipative relation for q as the Fourier equation and the 
dissipative relations for a as the definition of a Newtonian fluid; substituting the latter 
into conservation of momentum (3.8) yields 

p(v + D$) = -Dp + /j D 2 v + (A + lno)D(D ■ v) (4.16) 

which is the Navier-Stokes equation. These are the classical equations describing heat 
conduction and viscosity respectively, with kq being the thermal conductivity, /j,q the shear 
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viscosity and Ao the bulk viscosity. Originally, these dissipative relations were discovered 
experimentally, with Ao = 0, but the above method shows that they may be loosely derived 
from the local second law. This seems to have been originally noticed by Eckart [12], who 
immediately generalised it to relativity [2]. 

Note that one of the dynamical equations is the local Gibbs equation (4.7), even 
though this was originally concerned only with thermostatics. Thus it has been implicitly 
assumed that, in the general thermodynamic case, the Gibbs equation, as well as the 
thermal and caloric equations of state, continue to take their thermostatic forms. Also, 
note that there is no need for a state-space Gibbs equation t?ds* = de* +pd(l/p); indeed 
it would be generally inconsistent with the material Gibbs equation (4.7). In this article, 
the material Gibbs equation and the relation ip = q/d will be called entropic relations, 
since they effectively define the entropic fields (s,(p), up to the usual ambiguity in s. 
In the kinetic theory of gases, the entropic fields are defined in terms of the molecular 
distribution, allowing the entropic relations to be derived [5]. 

In summary, the basic thermo hydro dynamic fields are (p,v,p,a,e,q,'d,s,<p), which 
constitutes twenty functions. There are five conservation equations, namely the first law 
and conservation of mass and momentum. (The first law may be regarded as conservation 
of heat, as will become clear in the relativistic theory). The system is completed by two 
equations of state for (p, e*), four components of the entropic relations for (s, <p) and nine 
dissipative relations for ((?, cr), the latter being chosen to imply the second law. 

The fields pair up neatly with the equations expressing them or their material deriva- 
tives in terms of the basic fields and their spatial derivatives, except for temperature and 
the first law; even this becomes a temperature propagation equation on substituting other 
relations. For instance, for an inviscid ideal gas at constant density p = po, the local first 
law reads 

c oP o$ = k D 2 $ (4.17) 

which is the diffusion equation. This describes infinitely fast propagation of temperature 
and is therefore physically implausible. There is a similar problem of infinitely fast propaga- 
tion implied by the Navier-Stokes equations [13]. Thus the standard local thermodynamics 
described so far, e.g. by de Groot & Mazur [5], needs to be modified. 

V. Extended thermodynamics 

The problem of infinite propagation speeds stems from the assumption that the entropy 
density s satisfies the Gibbs equation (4.7), originally intended to apply only to the thermo- 
static case. A generalisation due to Miiller [3] involves allowing the the entropy to depend 
on quadratic terms in the dissipative fields (q, a). This reduces to the foregoing theory in 
an approximation where the quadratic terms may be neglected, i.e. the new theory is a 
quadratic rather than linear approximation to non-equilibrium. Applying the second law as 
before yields modified dissipative relations containing relaxation terms that provide finite 
propagation speeds. Moreover, dissipative relations of exactly this type may be obtained 
in the kinetic theory of gases by taking sufficiently many moments of the Boltzmann equa- 
tion, as originally shown by Grad [8]. Miiller & Ruggeri [13] now describe this as extended 
thermodynamics of irreversible processes, reserving extended thermodynamics to describe 
a more systematic theory which assumes a further divergence-type balance law. 
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Since a distinction between entropy and thermostatic entropy is now required, it is 
convenient to use s* for the thermostatic specific entropy henceforth, with s being the 
entropy density as before. The equations displayed in the preceding subsection have been 
written with this distinction in mind, e.g. the second law (4.6) contains s and the Gibbs 
equation (4.7) contains s*. The thermostatic relation s = ps* is now modified to include 
quadratic terms in (q,cr), giving the general form 

s = ps* \p{b q q ■ q + b m w 2 + b q q : q). (5.1) 

Here the coefficients b are arbitrary apart from being non-negative, since entropy must 
reach a maximum in the thermostatic case: 

b {q ,^ } > 0. (5.2) 

Similarly, the entropy flux ip may differ from the thermostatic entropy flux qjd by quadratic 
terms in (q,cr), giving the general form 

<f = | - k^wq - k^q ■ q. (5.3) 

Then the entropy production may be written explicitly as 



s + sD-v + D- (p = -q- ( + b q pq + k q D ■ q + k^Dw 



- w (^—^- + b^pzu + k m D ■ q 
_ ^ _ + hqp ^ + kq D q^j (5.4) 

where the coefficients (6, k) have been assumed constant for the sake of simplicity. So the 
second law suggests modified dissipative relations 

q = -K (D$ + b q p<d 2 q + k^ 2 D ■ q + k^Dw) (5.5a) 
w = -\ (D ■ v + b^p$w + k^-dD ■ q) (5.56) 
q = -2p (D ®v-\(D-v)h + b,p$q + k^(D <g> q - \{D ■ q)h)) (5.5c) 

which yield the entropy production with the same quadratic form as (4.15). The key point 
here is that the relaxation coefficients bf introduce the time derivatives of the dissipative 
fields / = (q, w,q) into the dissipative relations, giving equations for / + tff, where the 
relaxation timescales are 

t q = Kobqpi} 2 (5.6a) 
tzu = A 6 ro/ o# (5.66) 
t q = 2p b q p'&. (5.6c) 
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This leads to finite propagation speeds, at least for linear perturbations [13]. The temper- 
ature propagation equation for an inviscid ideal gas at constant density now becomes 

copo(d + t q 0) =k D 2 $ (5.7) 

for constant t q . This is the telegraph equation, with finite propagation speed a/ Ko/copot q . 

In summary, Miiller [3] obtained a causal theory of thermodynamics simply by allowing 
the entropic fields (s, cp) to depend on quadratic terms in the dissipative fields (q, a). There 
has been no tampering with the conservation laws, the second law, the equations of state 
or even the Gibbs equation, provided the latter is taken in thermostatic form, i.e. in 
terms of the thermostatic specific entropy s*. This is often explained misleadingly as a 
generalisation of the Gibbs equation; in reality, the thermostatic Gibbs equation is still 
assumed but the entropy density is no longer assumed to be purely thermostatic. One 
might continue to suppose that the thermostatic s* is a state function, but the entropy 
density s can no longer be a state function, even for an ideal gas, thereby toppling what 
many textbooks erect as a principle of thermodynamics. 

VI. General relativity 

General relativity is based on a four-dimensional manifold, space-time, with a symmetric 
bilinear form g, the metric. Taking units such that the Newtonian gravitational constant 
and the speed of light are unity, the Einstein equation is 

G = 8nT (6.1) 

where G is the Einstein tensor of g and T is the energy tensor of the matter, which will be 
taken in their contravariant forms. The energy tensor is more fully described as the energy- 
momentum-stress tensor, due to the physical interpretation of its various projections in a 
given frame as the energy density, momentum density or energy flux, and stress. 
The contracted Bianchi identity, a purely geometrical identity, reads 

V-G = (6.2) 

where V is the covariant derivative operator of g. Therefore 

V-T = (6.3) 

which expresses energy-momentum conservation. A material model is specified by giving 
the general form of T in terms of the material fields, supplemented by any required equa- 
tions other than energy-momentum conservation. This paper gives such a procedure as 
a definition of thermodynamic matter. Any such procedure requires at least a first law 
of thermodynamics. Physically, the first law is an energy conservation equation, but it 
differs from conservation of energy in the sense of a time-component of energy-momentum 
conservation (6.3). The difference turns out to be conservation of mass. In the Newtonian 
limit, conservation of energy yields conservation of mass to leading order and the first law 
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as the next correction. Thus relativistic thermodynamics requires a concept of mass. In a 
microscopic description, this would be just the mass of the particles. 

In relativity, mass is described by an energy-momentum vector J which will be called 
the material current. Its magnitude is the density p and its direction is the velocity vector 
u, assumed timelike: 

J = pu (6.4a) 
u-u b = -l. (6.46) 

Here the sign convention is that spatial metrics are positive definite and b denotes the 
covariant dual with respect to g (index lowering). Similarly, ft will denote the contravariant 
dual (index raising). These accents are included for book-keeping and are irrelevant for 
many purposes. In relativistic kinetic theory of gases [7], J is the first moment of the 
molecular distribution, with T being the second moment. Conservation of mass and energy- 
momentum may then be derived from the relativistic Boltzmann equation. 
Conservation of mass is expressed simply by 

V • J = 0. (6.5) 

Written explicitly, 

= V-J = p + pV-u (6.6) 

where the overdot henceforth denotes the covariant derivative along u, f = u • V/, which 
is the relativistic material derivative, or comoving derivative along the flow. This has a 
similar form to non-relativistic conservation of mass (3.7). 

Mass is a local form of energy in relativity, so it has an energy tensor Tm which will 
be called the material energy tensor. This takes the standard form 

T M = pu®u (6.7) 

so that the relation 

V- J = -u- (V-T M ) b (6.8) 

shows that conservation of mass is consistent with energy-momentum conservation for this 
type of matter, known as dust. 

VII. Heat 

The proposed theory of relativistic thermodynamics is based on four main physical prin- 
ciples. The first principle is that (i) heat is a local form of energy. In general relativity, 
heat is therefore described by an energy tensor Tjj which will be called the thermal energy 
tensor. Thermodynamic matter will be taken to mean a form of matter for which mass 
and heat are the only local forms of energy: 

T = T M + T H . (7.1) 
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Comparing with a microscopic description, Th is the effective energy tensor produced by 
the random motion of particles around the average flow given by u, as can be made precise 
in relativistic kinetic theory of gases [7] . 

The physical meaning of the various components of the thermal energy tensor are now 
determined: 

e = u-T b H -u (7.2) 

is the thermal energy density, 

q = -±{T H ■ u b ) (7.3) 

is the thermal flux and 

r = ±T H (7.4) 
is the thermal stress, where _L denotes projection by the spatial metric 

h~ x = g + u b <g> u b (7.5) 

orthogonal to the fluid flow. Note that h denotes the contravariant form and h~ x the 
covariant form. Therefore 

T H = eu <g> u + 2u <g> q + r. (7.6) 

As the notation indicates, thermal energy density, flux and stress will be identified respec- 
tively with the internal energy density, heat flux and stress of the non-relativistic theory. 
Therefore the energy tensor of the thermodynamic matter is 

T = (p + e)u <g> u + 2u <g> q + r. (7.7) 

This suffices to determine the quasi-local first law of relativistic thermodynamics in spheri- 
cal symmetry, as is shown in the accompanying paper [1] . Written in terms of the effective 
density g = p + e, the energy tensor has the same form as that of Eckart [2], who defined 
specific internal energy as g/p plus an undetermined constant. This is consistent, since 
specific thermal energy is e/p. 

The above principle implies that (e, q, r) are different aspects of heat. In particular, 
one may say that thermal energy density or heat density e represents heat at rest, and q 
heat in motion, relative to the average flow. This furnishes the common-sense explanation 
for such simple physical processes as the heating, insulation and cooling of a body, namely 
that heat flows in, remains at rest, then flows out. This intuitive explanation is often 
mocked in thermodynamics textbooks, accompanied by the claim that heat flux is the 
only form of heat. In relativity, thermal flux is just one projection of a thermal energy 
tensor. 

Similarly, what is traditionally called internal energy is just thermal energy, or simply 
heat, as the term is used in everyday language. Moreover, the ambiguous zero of internal 
energy is fixed for thermal energy. Recall that the original meaning of internal energy 
derived from the fact that the measurable energies, due to heat flux and work, were not 
conserved. According to conservation of energy, there had to be some other form of energy, 
internal to the system. The experiments of Joule established that this internal energy 
satisfied simple laws for an ideal gas. This stop-gap concept of internal energy differs 
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intrinsically from that of thermal energy, yet they turn out to agree. That internal energy 
is really thermal energy should anyway be clear even from non-relativistic kinetic theory 
of gases [5], where internal energy is defined as the average kinetic energy of the particles 
in the centre-of-momentum frame. 

The second physical principle is (ii) conservation of mass (6.5). As energy- momentum 
conservation (6.3) is automatic in general relativity, conservation of mass is equivalent by 
(6.8) to conservation of heat 

u-(V-T H f = 0. (7.8) 

Written explicitly, 

= - u ■ (V • T H ) b = e + eV ■ u + V • q + q ■ ii b + r : (V <g> v 9 ). (7.9) 

This reduces to the non-relativistic first law (4.5) in the Newtonian limit, since the term in 
the acceleration ii disappears, as shown in the next section. Thus principle (ii) is equivalent 
to the relativistic first law. This reveals that the first law is conservation of energy minus 
conservation of mass, cf. Eckart [2] . One may say that the first law is conservation of heat 
for the thermodynamic matter. Note that this could be modified in the presence of other 
matter, which would also contribute to the total energy tensor T. 

Several points perhaps need to be stressed. Firstly, although kinetic theory of gases 
has been mentioned for comparison, the above theory of heat is manifestly independent of 
such microscopic theories. Secondly, conservation of mass is required by the first law of 
thermodynamics. In the perspective of general relativity, this is the primary distinction 
between thermodynamic matter and any other energy tensor T to be used in the Einstein 
equation. Thirdly, conservation of mass involves a material current which determines a 
preferred time direction u, the material frame of Eckart. It is only in this frame that ther- 
mal flux — meaning the energy-flux component of the thermal energy tensor — agrees with 
what is normally called heat flux, even in the limit of non-relativistic thermodynamics. 
Finally, the basic conservation equations (6.3) and (6.5) are standard in relativistic ther- 
modynamics and date back to Eckart [2] . What seems to have been lacking is the physical 
interpretation in terms of mass and heat. In particular, conservation of heat should not 
be confused with the old caloric theory, as discussed further in the Conclusion. 

VIII. Newtonian limit 

The Newtonian limit may be described in terms of a vector A which plays the role of 
Newtonian time. The function $ defined by 

e 2<I> = -A-A b (8.1) 

turns out to reduce to the Newtonian gravitational potential. Similarly, the component of 
u orthogonal to A turns out to reduce to the Newtonian velocity: 

v = u + e~ 2$ (-u • A b )A 

Inverting, 

e*(-u — v) 

" \V\ 2 
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(8.2) 



(8.3) 



Taking units of length, factors of the speed c of light may be introduced by the formal 
replacements 



h i-> h (8.4a) 

(u, A, v) I— > c _1 ('u, A, u) (8.46) 

($,p)^c- 2 ($,p) (8.4c) 

(e,r) i-> c _4 (e,r) (8.4d) 

q i — c~ 5 q. (8.4e) 

These factors are determined by the desired physical interpretation of the fields in Newto- 
nian theory. This implicitly includes consequences like D i— > D for the covariant derivative 
of a -1 and / h- > c -1 / (where />—>/) for the material derivative f = u ■ V/. These 
replacements will be assumed in any subsequent equation involving c. 
For A to qualify as a candidate for Newtonian time, it should satisfy 

±L x h = 0(c~ 2 ) (8.5) 

where L\ is the Lie derivative along A. This condition reflects the absoluteness of Newto- 
nian space. Then 

±(X7 ®u b ) = D®v b + 0(c- 2 ) (8.6) 

and in particular 

V -u = D -v + 0(c- 2 ). (8.7) 
So conservation of mass (6.6) and heat (7.9) read 

p + pD -v = 0(c" 2 ) (8.8a) 
e + eD ■ v + D ■ q + r : (D <g> v b ) = 0(c~ 2 ) (8.86) 



which manifestly reduce to non-relativistic conservation of mass (3.7) and heat (4.5). 
Conservation of momentum may be written explicitly as 

= _L(V • T) = (p + e)±ii + _L(g • V)u + q(V • u) + ±q + i_(V • r). (8.9) 

Inserting factors of c as in (8.4), all terms except the first and last disappear in the Newto- 
nian limit. A longer argument, skipped here, confirms that this reduces to non-relativistic 
conservation of momentum (3.8). Thus all the conservation equations of non-relativistic 
thermodynamics have been generalised to general relativity, as conservation of energy- 
momentum and heat or mass. Similarly, it is straightforward to check the Newtonian 
limits of the second law, entropic relations and constitutive relations yet to be given. 

IX. Thermostatics and thermodynamics 

The next step, giving equations of state, requires a distinction between thermostatic and 
dissipative quantities, e.g. between thermostatic pressure and viscous pressure. Thus the 
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thermal energy tensor T H is divided into a thermostatic part T and a dissipative part T\ 
by 

T H = T + T 1 (9.1) 

with Ti = in the thermostatic case. The two are to be distinguished as in the non- 
relativistic theory by the type of equations determining them: To is to be given by equations 
of state and Ti by dissipative relations consistent with the second law, the latter being 
described in the next section. 

A fluid may be defined by T being isotropic, therefore taking the general form 

T = eu®u+ph. (9.2) 

Then p is the thermostatic pressure and e the thermostatic energy density. This isotropic 
form is appropriate for a fluid but not a solid; for instance, an elastic solid generally has 
anisotropic thermostatic stress, i.e. a general tensor replacing ph in To. 

The desired thermal and caloric equations of state may then be combined into an 
equation of state for To in terms of (p, #) . For instance, an ideal gas may be defined by 

T = (c u ®u + R h)p$ (9.3) 

which gives thermal and caloric equations of state with the same form as for the non- 
relativistic ideal gas (4.8). Generally, the equations determining T may be called the 
thermostatic relations. 

In the thermostatic case T± = 0, this recovers the relativistic model of a perfect fluid 
[14] , consisting of the energy tensor 

T M + T = (p + e)u <g> u + ph (9.4) 

together with equations of state for (e,p) as functions of p and conservation of mass or 
heat (7.9), which yields 

e* = pp/p 2 (9.5) 
in terms of the specific thermostatic energy 

e* = el p. (9.6) 

For an ideal gas, this gives i? = fop R °' c ° for a time-independent function f . Thus 

e = c f p 1+R " /co (9.7a) 
p = R foP 1+Rt)/co . (9.76) 

These equations of state specify a thermostatic ideal gas. Actually, it is more common to 
use the effective density g = p + e instead of p and e, forget mass conservation and give 
a single equation of state for p as a function of q. Indeed, it is common to ignore e and 
replace p with q in (9.7b), giving, for constant /o, the polytropic equation of state. This is 
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a reasonable approximation if e <C p, as in Newtonian theory, but the distinction between 
p and q is important in general. 

Returning to the general thermodynamic case, defining the dissipative stress or viscous 
stress 

a = T -ph (9.8) 

and the dissipative energy density 

r] = e-e (9.9) 
allows the dissipative energy tensor to be composed as 

T x = -qu <g> u + 2u <g> q + a (9.10) 

which is the general form of such a tensor. It will turn out that rj vanishes in the linear 
approximation for entropy, allowing consistency with standard non-relativistic thermody- 
namics [5]. Nevertheless, rj will be retained on the grounds that, as a matter of principle, 
thermal energy need not be purely thermostatic. The same distinction between thermo- 
static and dissipative energy density could be drawn in the non-relativistic case, leaving 
ambiguities to be fixed in equations where rj might occur; these were fixed in equations 
displayed in Section IV by using e and e* with appropriate foresight. 
Exhibiting all terms in the energy tensor, 

T = T M + T + Ti = (p + e + rj)u u + 2u <g> q + ph + a. (9.11) 

This shows that T incorporates all the basic thermohydrodynamic fields of the non- 
relativistic theory described previously, apart from temperature $ and entropy (s,cp), but 
adding dissipative energy density rj. 

X. Entropy 

The third physical principle is that (iii) entropy is a local current. In relativity, it is 
therefore described by an entropy current vector \& whose components are entropy density 

s = -u b -y (10.1) 

and entropy flux 

ip = _Ltf. (10.2) 

As their names indicate, entropy density and entropy flux will be identified with their 
non-relativistic versions. Therefore 

$> = su + ^. (10.3) 

The fourth physical principle is (iv) non-destruction of entropy. In relativity, this is there- 
fore expressed by 

V-*>0. (10.4) 
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Explicitly, 



s + sV-u + V-(p>0 



(10.5) 



which reduces to the non-relativistic second law (4.6) in the Newtonian limit. Thus prin- 
ciple (iv) is the relativistic second law. These two principles are standard in relativistic 
thermodynamics. Again it seems to be Eckart [2] who first gave a version of (10.5), as- 
suming <p = q/'d. 

The final stage consists of introducing entropic relations for \1> in terms of the other 
fields, then finding dissipative relations for T\ which imply the second law. Successive 
approximations for \1/ will be taken: \1/ = \l/o for the thermostatic case, ^ = ^0 + ^1 
for the linear approximation and \& = ^0 + ^1 + ^2 for the quadratic approximation for 
entropy. 

The thermostatic entropy current W may be assumed to satisfy 

= (10.6a) 
V-* = -^-(V-T ) b . (10.66) 

The first equation states that the thermostatic entropy current is comoving with the ma- 
terial. Therefore \I/o may be written in terms of the specific thermostatic entropy 

s* = -u b -#o/p- (10.7) 

For a fluid (9.2), the second equation may be written as 

#V-#o = e + eV-u + pV-u (10.8) 

or 

$s*=e*-pp/p 2 . (10.9) 

This form was given by Eckart [2] and is identical to that of the non-relativistic local Gibbs 
equation (4.7).* This is usually regarded as defining thermostatic temperature. For an 
ideal gas (9.3), it integrates to 

s* = c log #-#0 log p (10.10) 

which has the same form as in the non-relativistic case (4.9). 

* Most references on relativistic thermodynamics give a different Gibbs equation. Some- 
times the material derivative is replaced with the covariant derivative in the Gibbs equation 
[15,16], giving an equation whose spatial components are quite unnecessary and generally 
overdetermine the fields. This equation has even been described as the first law. In other 
cases, the material derivative is replaced with a state-space differential in the Gibbs equa- 
tion [4,6], which generally does not imply the material Gibbs equation (10.9). It has 
even been claimed that the thermostatic (^o 7 Tq) should be replaced with (\I/,T) in this 
state-space equation [4]. 
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The relations (10.6) therefore constitute a generalised Gibbs equation, relativistically 
unified in terms of \&o an d T . This relativistic Gibbs equation has the desired property 
that entropy is conserved in the thermostatic case: 

Ti = V • = P«* = (10.11) 

as follows from conservation of heat (7.8). In other words, thermostatics requires equality 
in the second law. Indeed, one might propose a stronger version of the second law, to the 
effect that non-thermostatic processes necessarily produce entropy. One could then define 
thermostatics, or equilibrium thermodynamics, or reversible thermodynamics, by equality 
in the second law (10.4). 

XI. Linear approximation 

The linear correction to entropy current may be assumed as 

V^-j-Tt (11.1) 

which may be regarded as defining temperature $ away from equilibrium. Then the entropy 
flux = q/fi has the same form as in non-relativistic thermodynamics. The component 
of along the flow, 77/$, is not so determined, except that the above linear form is the 
simplest. However, it turns out to yield a useful cancellation in the entropy production, 
which, using conservation of heat (7.8), becomes simply 

V • (#0 + *i) = -Ti : (V <g> j^j . (11.2) 

Thus the entropy production is a contraction of the dissipative energy tensor T\ with 
the gradient of the inverse temperature vector u/fi, thereby unifying what are sometimes 
called, respectively, the thermodynamic fluxes and thermodynamic forces [5]. Writing the 
components explicitly, 

y.(» + »0 = - «-<™ + ^> - g! <7 tt> > -:g. (11.3) 

So the second law (10.4) suggests the dissipative relations 

q = -«o-L(V tt i? + &u) (11.4a) 

zu = -X V-u (11.46) 

q = -2^o^(V H ®u- l(V-u)h) (11.4c) 

rj = -z/ tf (11.4d) 

with 

kq > A > /u > i/ > 0. (11.5) 



19 



If these coefficients are non-zero, the entropy production is given explicitly by 



The dissipative relation for q is a relativistic modification of the Fourier equation, including 
a term in the acceleration u which disappears in the Newtonian limit. The dissipative 
relations for a are relativistic versions of the viscosity relations for a Newtonian fluid, 
yielding a relativistic Navier-Stokes equation. With w = 0, these were originally derived 
by Eckart [2]. The last dissipative relation for r\ appears to be new and modifies the 
temperature propagation equation: taking the example of an inviscid ideal gas at constant 
density with vanishing acceleration, the first law (7.9) becomes 

PqCqV - = k D 2 $. (11.7) 

Thus the parabolic (diffusion) equation for vq = becomes a hyperbolic (telegraph) equa- 
tion for vq < 0, as desired, but an elliptic equation for v > 0. Since ellipticity is even 
worse than parabolicity on causal grounds, this leads to v Q = and therefore r\ = as 
promised. Thus the linear approximation for entropy recovers and generalises the Eckart 
theory [2] , which concerned the case of a fluid without bulk viscosity. 

XII. Quadratic approximation 

For the quadratic correction ^ 2 to the entropy current, the most general form which is 
quadratic in the dissipative energy tensor T\ is 

^2 = • 0° + b^w 2 + : ^ + b v rj 2 )pu - k m wq - k^q ■ q b - k v rjq (12.1) 

where 

b {g ,^M > ( 12 - 2 ) 

so that entropy is maximised in equilibrium. The coefficients b lead to relaxation effects 
as in the non-relativistic case. Explicitly, the entropy production is 

_|_ flyb 

^2 + b g pq b + K(V-s) b + k^Vw + k v Vr) 

- & {—J- + b^pw + k m V ■ q 

-vl^ + b^PV + k^-q) (12.3) 
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where the coefficients (6, k) have been assumed constant for the sake of simplicity, the 
generalisation being straightforward. The entropy production takes the same quadratic 
form (11.6) as in the linear approximation if the dissipative relations are taken as 

q = -Ko-L(V tt # + $ii + b q p$ 2 q + k^ 2 V ■ q + k^V^w + k^V^rf) (12.4a) 

w = -A (V- u + b m p$vj + k^dV-q) (12.46) 

q = -2/j ±(V tt <g> u - |(V • u)h + b,p$q + fc s tf(V* <g> q - \h : (V <g> o b )/i)) (12.4c) 

7/ = -z/ (0 + ^ptf 2 r) + k v $ 2 V ■ q) . (12 Ad) 

This completes the system of thermodynamic field equations for specified coefficients (b, k). 
In the case r\ = 0, these dissipative relations reduce to those of Israel & Stewart [4], also 
given by Jou et al. [17]. If the coefficients (6, k) are not constant, additional derivatives 
appear, as noted by Hiscock & Lindblom [18]. 

Relativistic kinetic theory of gases [7] confirms the existence of dissipative energy 
density rj; indeed it emerged unrecognised in the non-relativistic theory of Grad [8] as 
essentially the 14th moment of the molecular distribution, the previous 13 moments being 
(p, v, e, q, q), with being given by equations of state and w vanishing. For instance, 

for an monatomic (4.11) ideal (9.3) gas with w = 0, setting 

&-— k-— &-— k--^- b--^- (12 5) 

and 

hp 2 2p 15p 3 

K0= 2B^ ^ = W^p VQ= mpH (12 - 6) 
where So is a constant, the dissipative relations in the non-relativistic limit become 

= q + B pq + \l$p&D<& + Rq-OD ■ q + ^Drj (12.7a) 
= q + lB pq + 2R p$(D <g> v - \{D ■ v)h) + §(£> <g> q - \{D ■ q)h) (12.76) 
= 77 + B prj + Rq-OD ■ q (12.7c) 

where the speed of light has been introduced as in (8.4), with Rq i— > c 2 Rq, Bq i— > cBq, 
$ i— > c -4 ^ and rj i— > c~ 6 r/. Together with the non-relativistic conservation laws, these are 
the 14-moment equations (1.42) of Miiller & Ruggeri [13], their 14th moment being 877, 
with the factor determined in relativistic kinetic theory [7] . This confirms the existence of 
dissipative energy density and thereby topples another pillar of textbook thermodynamics, 
that internal energy necessarily be a state function. 

XIII. Integral quantities and laws 

There exist integral forms of the second law and conservation of mass, because the local 
forms are current (semi-) conservation laws, allowing the Gauss divergence theorem to be 
applied. In particular, the mass of a region E of a spatial hypersurface may be defined by 



M(E) = - I J (13.1) 
is 
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where the volume form and contraction with unit normal are implicit. Consider a space- 
time region O bounded by E, a later E' and a hypersurface region E generated by flowlines 
of u. Applying the Gauss theorem, 

M(E')-M(E)= f V-J = (13.2) 

which is the integral form of conservation of mass. Similarly, the entropy of E may be 
defined by 

5(E) = -J * (13.3) 

and the entropy supply through E by 

So(E) = - / (13.4) 
J± 

Then the Gauss theorem yields 

S(E') - 5(E) - S {±) = [ V • * > (13.5) 

which is the integral second law. Thus there is an integral second law of thermodynamics 
even in general relativity. An integral first law may seem awkward in similar generality, 
since the local first law (7.8) is not a current conservation law, but in spherical symmetry 
a simple integral (or quasi-local) first law exists, as shown in the accompanying paper [1]. 
The integral laws become particularly simple if the vorticity vanishes, 

±(VAw b )=0 (13.6) 

where A denotes the antisymmetric tensor product. Then there exist spatial hypersurfaces 
orthogonal to u which may be used to define preferred volume integrals. For instance, the 
mass is 

M= I *p (13.7) 



where E henceforth denotes a region of a hypersurface orthogonal to u, with volume form 
*. The induced metric of E is just which defines *. One may assume that E has finite 
volume, or infinite volume with all relevant integrals existing. 

Concerning time derivatives: in the relativistic case, M is generally meaningless, since 
proper time along u generally does not label the hypersurfaces orthogonal to u. Labelling 
these hypersurfaces by t, choosing vanishing shift vector and denoting Af = df/dt for 
functions /, the relationship 

/ = e"*A/ (13.8) 

defines the function 0, up to relabellings of t. Then AM is meaningful and the Gauss 
theorem yields 

V- J = I *p- I *p (13.9) 
n Jn' Jt, 
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which differentiates to 

! *e V-J = AM. (13.10) 
Jt 

Therefore local conservation of mass (7.10) implies the integral form 

AM = (13.11) 

which has the same form as classical conservation of mass (2.8). 
Similarly, the entropy is 

S = I *s (13.12) 
Jt 

and the entropy supply through E is 

S = - [ *-<p (13.13) 
Jt 



where * is the volume form of E times the unit outward normal covector tangent to E. 
Then 

'as 



>o = - / • e"V (13.14) 

JdT 



where the surface <9E is the intersection of E and E, with area form and unit normal implicit 
in §. This has the same form as the classical definition of entropy supply (4.4), apart from 
the use of A and the corresponding factor e^, which tends to one in the Newtonian limit. 
Applying the Gauss theorem again, 



which differentiates to 



V-^=/ *s- *s+ *-<p (13.15) 
n Jt,' Jt, Jt, 



[ * e ^V- ^ = AS- AS . (13.16) 



Therefore V • * is the entropy density production rate and the local second law (10.4) 
integrates to 

AS > AS (13.17) 

which has the same form as the classical second law (2.5). 
Similarly, the heat supply through E may be defined by 



Q = - *-q (13.18) 
Jt 

which yields 

AQ = -(f -e^q (13.19) 

JdT, 
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again with the same form as the classical definition of heat supply (4.3). Finally, the 
thermal energy may be defined by 



explain the notation adopted earlier: the material and thermal energy tensors Tm and T# 
are the energy tensors of mass and heat respectively. 

Thermal energy H is what thermodynamicists call internal energy and what ordinary 
people call heat. In particular, this defines the heat of a body or volume of fluid. Ther- 
modynamics textbooks often claim that there is no such thing as the heat of a body, in 
bizarre conflict with common experience. This seems to involve the belief that heat flux 
is the only form of heat, i.e. what was really meant is that there is no conserved quantity 
corresponding to the flux of heat, as supposed by the old caloric theory. A nice analogy 
[11], intended to support the above claim but actually revealing the opposite, is that there 
is no such thing as the amount of rain in a lake. True, because rain is not the only form 
of water, but the amount of water in the lake does make sense. 

XIV. Conclusion 

A general macroscopic theory of generally relativistic thermodynamics has been developed, 
based on clear principles and sufficient to include relaxation effects. There is a full set of 
dynamical equations, easily applied in practice, with no meaningless derivatives or other 
common thermostatic confusions. The basic fields are the temperature the material, 
thermostatic and dissipative energy tensors (Tm, Tq, Ti) and the entropy current These 
fields may be paired respectively with the thermodynamic field equations: conservation 
of mass (6.5) or heat (7.8), conservation of energy-momentum (6.3), the thermostatic 
relations for Tq, e.g. the equations of state (9.3) for an ideal gas, the dissipative relations 
for Ti and the entropic relations for \1>: thermostatic (10.6), linear (11.1) and quadratic 
(12.1) approximations. The dissipative relations (11.4) or (12.4) are chosen to yield the 
second law (10.4), subject to relativistic causality. Having ensured these conditions, the 
dynamical system may be reduced by eliminating \1/ and the entropic relations, since W 
does not appear in the remaining equations. Similarly, the thermostatic relations may be 
used to eliminate further fields, numbering two in the case of a fluid. 

Four physical principles have been emphasised, which may be regarded as minimal 
requirements for a definition of thermodynamic matter. In principle this includes solids, 
necessarily imperfect in relativity, though the form of the thermostatic relations has been 
specified only for fluids. Three of the principles are standard. However, the first principle, 
that heat is a local form of energy, seems not to have been clearly stated before. This 
principle unifies various thermodynamic quantities in an essentially relativistic way as the 




(13.20) 



The forms 




(13.21a) 



(13.216) 
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energy, flux and stress components of a thermal energy tensor. This unification is analogous 
to that of electricity and magnetism in relativistic electromagnetism. Moreover, this agrees 
with non-relativistic thermodynamics in the physical interpretation of thermal flux only 
in the material frame. The principle thereby resolves the controversy of Eckart [2] versus 
Landau-Lifshitz [6] frame. 

A further consequence of this formulation of basic principles is that it is unnatural 
to insist that thermal energy be purely thermostatic, despite the traditional dogma that 
internal energy be a state function. The similarly traditional dogma that entropy be a 
state function had already disappeared in Miiller's extended thermodynamics [3]. Thus we 
have an apparently new thermo dynamical quantity, dissipative energy, which nevertheless 
occurred unrecognised in the kinetic theory of gases [7]. These former assumptions may 
now be derived in the linear approximation to non-equilibrium, but need not hold in 
general, e.g. in the quadratic approximation. 

One artificial feature of this theory, shared with the original Miiller [3] and Israel- 
Stewart [4] theories, is that the quadratic correction to entropy current involves certain 
coefficients which are left unspecified. However, these coefficients can be determined for 
an ideal gas by kinetic theory of gases [13]. Kinetic theory also automatically implies 
conservation of mass and energy-momentum, and higher moments can be used to obtain 
a hierarchy of divergence equations for which there is a standard mathematical theory of 
local existence, uniqueness and causality [19]. However, this type of relativistic extended 
thermodynamics has so far been worked out under the assumption of vanishing dissipative 
energy density [13]. The main physical problem with this type of theory is that it is not 
known how to directly measure the higher moments. 

It seems necessary to stress that conservation of heat for thermodynamic matter has 
been derived from conservation of mass and energy-momentum, and definitely differs from 
the old caloric theory of Carnot [20]. A relativistic version of the caloric theory would be 
to describe heat by an energy-momentum vector, the caloric current Jh = su + q, and 
conservation of heat by V • Jh = 0, as for conservation of mass. This would integrate 
by (13.18) and (13.20) to AH = AQ. Carnot's theory of heat can be interpreted as in- 
volving this equation and the claim that the heat if is a state function, which contradicts 
experiment. Unfortunately the discrepancy between this and the classical first law (2.1) 
has led to a historical over-reaction against the intuitive notion of heat, specifically the 
widespread insistence that the internal energy of Clausius differs from the heat of Carnot, 
whereas Carnot's mistake was instead in the formulation of conservation of heat. The 
relativistic theory reveals that the caloric theory failed only because the mathematical na- 
ture of energy — an energy-momentum-stress tensor Th rather than an energy-momentum 
vector Jh — was not understood in pre-relativistic physics. In particular, material stress 
appears in the form of work in the first law, because it is also a form of heat. If the 
thermodynamic matter is not isolated, heat and other forms of energy may be exchanged. 

The relativistic theory similarly resolves the ancient philosophical question of whether 
heat is a substance. Put simply, both mass and heat are forms of energy, but with different 
natures: mass is described by a current J which satisfies a vectorial conservation law (6.5), 
conservation of mass, whereas heat is a general form of energy Th satisfying a tensorial 
conservation law (7.8), the first law of thermodynamics. 
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